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Abstract 

In this paper we briefly review the main aspects of a recent pro- 
posal to simulate semiclassical corrections to classical dynamics by 
suitable classical stochastic fluctuations, and we apply it to the spe- 
cific instance of charged beams in particle accelerators. The resulting 
picture is that the collective beam dynamics, at the leading semiclassi- 
cal order in Planck constant can be described by a particular diffusion 
process, the Nelson process, that is time-reversal invariant. Its dif- 
fusion coefficient ^/N\ c represents a semiclassical unit of emittance 
(here N is the number of particles in the beam, and A c is the Comp- 
ton wavelength). The stochastic dynamics of the Nelson type can be 
easily recast in the form of a Schroedinger equation, with the semi- 
classical unit of emittance replacing Planck constant. Therefore we 
provide a physical foundation to the several quantum-like models of 
beam dynamics proposed in recent years. We also briefly touch upon 
applications of the Nelson and Schroedinger formalisms to incorporate 
the description of coherent collective effects. 



1 To appear in the Proceedings of the International Workshop on "Quantum Aspects 
of Beam Dynamics", held in Stanford, 4-9 January 1998. 



1 



1 Introduction 



The dynamical evolution of beams in particles accelerators is gov- 
erned by external electromagnetic forces and by the interaction of the 
beam particles among themselves and with the environment. Charged 
beams are therefore higly nonlinear dynamical systems, and most of 
the studies on colliding beams rely either on classical phenomena such 
as nonlinear resonances, or they are concerned with isolated sources 
of unstable behaviors as building blocks of more complicated chaotic 
instabilities. 

This line of inquiry has produced a general qualitative picture of 
dynamical processes in particle accelerators at the classical level. How- 
ever, the coherent oscillations of the beam density and profile require, 
to be explained, some mechanism of local correlation and loss of statis- 
tical independence. This fundamental observation points towards the 
need to take into account all the interactions as a whole. Moreover, 
the overall interactions between charged particles and machine ele- 
ments are really nonclassical in the sense that of the many sources of 
noise that are present, almost all are mediated by fundamental quan- 
tum processes of emission and absorbtion of photons. Therefore the 
equations describing these processes must be, in principle, quantum. 

Starting from the above considerations, two different approaches 
to the classical collective dynamics of charged beams have been de- 
veloped, one relying on the Fokker-Planck equation for the beam 
density, another based on a mathematical coarse graining of Vlasov 
equation leading to a quantum-like Schroedinger equation, with a 
thermal unit of emittance playing the role of Planck constant || . 

The study of statistical effects on the dynamics of electron (positron) 
colliding beams by the Fokker-Planck equation has led to several in- 
teresting results, and has become an established reference in treating 
the sources of noise and dissipation in particle accelerators by standard 
classical probabilistic techniques (3| . 

Concerning the relevance of the quantum-like approach, at this 
stage we only want to point out that some recent experiments on 
confined classical systems subject to particular phase-space boundary 
conditions seem to to be well explained by a quantum-like (Schroedinger 
equation) formalism Q. 

In any case, both approaches do not take into account quantum 
corrections, while in principle these effects should be relevant, expe- 
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cially in fixing fundamental lower limits to beam emittance. In this 
report we give a short summary of a recently proposed model for the 
description of collective beam dynamics in the semiclassical regime. 
This new approach relies on the idea of simulating semiclassical cor- 
rections to classical dynamics by suitable classical stochastic fluctu- 
ations with long range coherent correlations, whose scale is ruled by 
Planck constant ||. 

The fluctuative hypothesis has been introduced by simple stability 
criteria, and it has been semiquantitatively tested for many stable sys- 
tems, including beams. The virtue of the proposed semiclassical model 
is twofold: on the one hand it can be formulated both in a probabilis- 
tic (Fokker-Planck) fashion and in a quantum-like (Schroedinger) set- 
ting. It thus bridges the formal gap between the two approaches. At 
the same time it goes further by describing collective effects beyond 
the classical regime due to the semiclassical quantum corrections. 

In particular, implementing the fluctuative hypothesis qualitatively 
by simple dimensional analysis, we derive a formula for the phase- 
space unit of emittance that connects in a nontrivial way the number 
of particles in the beam with Plank constant. 

The fluctuative scheme is then implemented quantitatively by in- 
troducing a random kinematics in the form of a diffusion process in 
configuration space for a generic representative of the beam (collective 
degree of freedom). 

We are interested in the description of the stability regime, when 
thermal dissipative effects are balanced on average by the RF energy 
pumping, and the overall dynamics is conservative and time-reversal 
invariant in the mean. Therefore, we model the random kinematics 
with a particular class of diffusion processes, the Nelson diffusions, 
that are nondissipative and time-reversal invariant (We will briefly 
comment at the end of the last section on the extension of the present 
scheme to include the treatment of dissipative effects) . 

The diffusion process describes the effective motion at the meso- 
scopic level (interplay of thermal equilibrium, classical mechanical sta- 
bility, and fundamental quantum noise) and therefore the diffusion 
coefficient is set to be the semiclassical unit of emittance provided 
by qualitative dimensional analysis. In other words, we simulate the 
quantum corrections to classical deterministic motion (at leading or- 
der in Planck constant) with a suitably defined random kinematics 
replacing the classical deterministic trajectories. 
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Finally, the dynamical equations are derived via the variational 
principle of classical dynamics, with the only crucial difference that the 
kinematical rules and the dynamic quantities, such as the Action and 
the Lagrangian, are now random. The stochastic variational principle 
leads to a pair of coupled equations for the beam density and the 
beam center current velocity, describing the dynamics of beam density 
oscillations. It is an effective description in the stability regime. 

The stochastic variational principle for Nelson diffusions (with dif- 
fusion coefficient equal to Planck constant) is a well developed mathe- 
matical tool that has originally been introduced to provide a stochastic 
formulation of quantum mechanics. Therefore, apart from the differ- 
ent objects involved (beam spatial density versus Born probability 
density; Planck constant versus emittance), the dynamical equations 
of our model formally reproduce the equations of the Madelung fluid 
(hydrodynamic) representation of quantum mechanics. In this sense, 
the present scheme allows for a quantum-like formulation equivalent 
to the probabilistic one. 

At the end of the last section we will briefly discuss how the hy- 
drodynamic formulation of the equations for the collective stochastic 
dynamics can be used to control the beams, for instance by selecting 
the form of the external potential needed to obtain coherent oscilla- 
tions of the beam density. 

2 Simulation of semiclassical effects by 
classical fluctuations 

Let us consider a physical system subject to a classical force law of 
modulus F(r) that is attractive and confining at least for some finite 
space region with a characteristic linear dimension R. Given iV ele- 
mentary granular constituents of the system, each of mass m, let v 
denote their characteristic velocity, and r their characteristic unit of 
time. 

A characteristic unit of action per particle is then defined as 

a = mv 2 T. (1) 

If the system has to be stable and confined, one must impose that 
the characteristic potential energy of each particle be on average equal 
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to its characteristic kinetic energy (virial theorem): 



C mv\ (2) 

where C is the work performed by the system on a single constituent. 
On the other hand, if the system extends on the characteristic length 
scale R, 

C NF(R)R. (3) 

By equations (2) and (3) we can express the characteristic velocity v 
as 

Inf(r)r 

v ~\ • 4 

Introducing the global time scale T associated to the system, we also 
have v = R/T. Replacing this expression and equation (4) for each 
power of v in equation (1), we obtain the following expression for the 
action per particle: 



mF(R)R 3 / 2 N 1 / 2 ^. (5) 

Mechanical stability requires that the action per particle must not 
depend on N, while on the other hand, the microscopic unit of time 
r must obviously depend on N and on the system's global time scale 
T. Therefore we must impose 

(6) 



N 

Inserting equation (6) into equation (5) we obtain the unit of action 
per particle as a explicit expression in terms of the constituent's mass, 
the system's linear dimension R, and the classical force calculated in 
R: 

a^m 1/2 R 3/2 ^F(R). (7) 

The scaling relation (6) can be also interpreted as a fluctuative hy- 
pothesis connecting the time scale of a microscopic stochastic motion 
with the classical time scale of the global system. In fact, equation 
(6) was first postulated by F. Calogero in his attempt to prove that 
quantum mechanics might be interpreted as a tiny chaotic component 
of the individual particles' motion in a gravitationally interacting Uni- 
verse 101. 
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In our scheme, rather than being a postulated consequence of clas- 
sical gravitational chaos, the fluctuative hipothesis of Calogero derives 
from a condition of mechanical stability. Since the stability conditions 
and the virial theorem apply to any classically stable and confined sys- 
tem, even with a small number of degrees of freedom, our derivation of 
equations (6) and (7) is universal as it applies to any interactions, not 
only gravity, and to systems composed by any number of constituents, 
not necessarily large, and not necessarily classically chaotic. 

We have verified that for any stable aggregate, plugging in equation 
(7) the pertaining interaction F, individual constituents' mass m and 
aggregate's linear dimension R, one has that the unit of action per 
particle a is always equal, in order of magnitude, to Planck action 
constant h. 

Our interpretation of this remarkable result is then that the fluc- 
tuative relation (6) and the associated formula (7) for the Planck 
quantum of action simulate (reformulate) in a classical probabilistic 
language the Bohr-Sommerfeld quantization condition. They provide 
a classical description of quantum corrections to classical phase-space 
dynamics at the leading semiclassical order h. 

We here briefly derive the result for the case of interest of a sta- 
ble bunch of charged particles in a particle accelerator. We consider 
a single electron (proton), in the reference frame comoving with the 
bunch. Confinement and stability of the bunch arise from the many 
complicated interactions among its constituents and between the same 
constituents and the external magnetic and RF fields. The net effect 
can be, in first approximation, schematized by saying that the single 
electron (proton) experiences an effective harmonic force, the typi- 
cal phenomenological law of force for beams when higher anharmonic 
contributions can be neglected: F(R) = KR, where K is the effective 
phenomenological elastic constant. We then have for beams: 

a * m^RiR 1 / 2 . (8) 

Let us consider for instance the transverse oscillations for protons 
at Hera: in this case we have K = W~ 12 Nm -1 , the linear transverse 
dimension of the bunch R = 10~ 7 m, and the proton mass. For elec- 
trons in linear colliders we have instead K = \Q~ ll Nm~ l , R = lCT 7 m, 
and the electron mass. In both cases, from equation (8) we have that 
the unit of action per particle a, ruling the coherence and stability of 
the bunch, is in both cases h, up to at most one order of magnitude. 
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All other instances of charged bunches considered lead to the same 
result, yielding our first important conclusion: the stability of charged 
beams is ruled by quantum effects on a mesoscopic scale. Moreover, at 
the semiclassical level, such quantum aspects can be described in terms 
of suitable classical fluctuations that mimick (simulate) the weak but 
unavoidable presence of fundamental quantum noise. 

The parameter that rules the stability of the system at the meso- 
scopic scale is however not directly h, but in the case of charged beams 
some characteristic unit of emittance. This is a scale of action, or of 
length when divided by the Compton wavelength, that measures the 
spread of the bunch in phase space, or, equivalently, in real space. 

This notion is very useful in the regime of stability and of thermal 
equilibrium that we explicitely consider. In this case the emittance can 
be expressed as a unit of equivalent thermal action. To introduce a 
characteristic unit of emittance in our fluctuative semiclassical scheme 
we then proceed as follows: the time scale of quantum fluctuations is 
defined as the ratio between h and a suitable energy describing the 
equilibrium state of the given system. This leads naturally to iden- 
tify this energy with the equivalent thermal energy fc^T, with ks the 
Boltzmann constant and T the equivalent temperature. On the other 
hand, in our scheme such time scale coincides with the fluctuative time 
r; we therefore have: 

h , , 

Using relation (6) we obtain the equivalent thermal unit of action 

ksTT = h\fN . (10) 

Introducing the Compton wavelength A c = h/mc and dividing by it 
both sides of equation (10) we finally obtain the characteristic unit of 
emittance £: 

£^X C ^/N. (11) 

Equation (11) connects in a nontrivial way the number of par- 
ticles in a given charged beam and the Compton wavelength. The 
square root of TV appears as a semiclassical "memory" of quantum 
interference. The relation (11) seems to point out the existence of a 
mesoscopic lower bound on the emittance some orders of magnitude 
above the quantum limit given by the Compton wavelength. More- 
over, Equation (11) yields the correct order of magnitude in for the 
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emittance in typical accelerators: for instance, with N = 10 -j-10 3 
one has £ = 10 _6 m in excellent agreement with the lowest emittances 
that are at the moment experimentally attainable. 

Actually, limits and requirements on beam existence, luminosity 
and statistics do not allow for beams with a number of particles ap- 
preciably lower than N = 10 10 -j- 10 11 . Thus the estimate (11) re- 
ally provides an a priori lower bound, as it implies that the emit- 
tance cannot be reduced appreciably below the mesoscopic thresholds 
£ = 10 5 -T- 10 6 A C , well above the Compton wavelength limit and only 
one or two orders of magnitude below the current experimental limits. 
It seems also unlikely that further quantum corrections beyond the 
leading semiclassical order could somehow contribute in lowering the 
mesoscopic bound (11) function of N. 

3 Stochastic collective dynamics in the 
stability regime 

The previous discussion can be made more quantitative by observing 
that the fluctuative relation (6) can be be recast with a little bit of 
work in the alternative form 

I ~ r 2 / 3 , (12) 

where I is a characteristic mean free path per particle. The detailed 
derivation of relation (12) from equation (6) is reported elsewhere 
||. Relation (12) indicates that the classical flcutuative simulation 
of semiclassical corrections really implies a fractal space-time relation 
in the mean, with a Kepler exponent associated to stable, confined 
and coherent dynamical systems, for instance charged beams in the 
stability regime. 

We therefore model the spatially coherent fluctuations (6) and (12) 
by a random kinematics performed by some collective degree of free- 
dom q(t) representative of the beam. The most universal continuous 
random kinematics that we can choose is a diffusion process in real 
or configuration space. In this way the random kinematics provides 
an effective description of the space-time variations of the particle 
beam density p(x, t) as it coincides with the probability density of the 
diffusion process performed by q(t). 



8 



Since it measures a collective effect at the mesoscopic scale, the 
diffusion coefficient must be related to the equilibrium parameter in 
the stability regime, that is to the characteristic semiclassical unit of 
emittance (11) rather than to the Plank action constant. 

Then, in suitable units, the basic stochastic kinematical relation is 
the Ito stochastic differential equation 

dq(t) = v+(q, t)dt + £ 1/2 dw , (13) 

where v+ is the deterministic drift, the square root of the characteristic 
emittance (11) is the diffusion coefficient, and dw is the time increment 
of the standard (5-correlated Wiener noise. 

We are concerned with the regime of stability of the beam oscil- 
lation dynamics, both since it is the relevant regime in the physics of 
accelerators and because the beam can be considered quasistationary 
during it, until, eventually, space charge effects become dominant and 
the beam is lost. In such stationary regime the energy lost by photonic 
emissions is regained in the RF cavities, and on average the dynamics 
is still time-reversal invariant. We can therefore still define a classical 
Lagrangian L(q, q) for the system, however with the classical deter- 
ministic kinematics replaced by the random diffusive kinematics (13). 

The equations of dynamics can then be deduced from the classical 
Lagrangian by simply modifying the variational principles of classical 
mechanics into stochastic variational principles. In fact, the mathe- 
matical techniques of stochastic variational principles have been de- 
veloped and applied to obtain Nelson stochastic mechanics, an in- 
dependet stochastic reformulation of quantum mechanics in terms of 
time-reversal invariant Markov diffusion processes with diffusion coef- 
ficient h j/J . In the context of Nelson stochastic mechanics one derives 
Schroedinger equation in the form of the Madelung coupled hydrody- 
namic equations for the probability density and the probability current 

In the present mesoscopic context the analysis is quite similar 
to that of Nelson stochastic mechanics, yielding again two coupled 
nonlinear hydrodynamic equations, however, with the emittance (11) 
replacing Planck constant in the diffusion coefficient, the real space 
bunch density replacing the quantum mechanical probability density, 
and the bunch center velocity replacing the quantum mechanical prob- 
ability current. 
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We now briefly sketch the derivation of the dynamical equations. 
The detailed analysis may be found elsewhere f8|. Given the stochastic 
differential equation (13) for the diffusion process q(t) in d = 3 space 
dimensions, one introduces the classical Lagrangian 

L(q,q) = l/2mq 2 -V(q). (14) 

For the generic trial diffusion q{t) one has, respectively, the prob- 
ability density p(x,t), the forward drift v + (x,t) and the backward 
drift V-(x,t). It is then useful to define two new variables, v(x,t) 
and u(x,t), respectively the current velocity and the osmotic velocity, 
defined as: 

v++v- v+-v- Vp , , 

U = ^2— 5 U= ^2~ =£ —- (15) 

The mean classical action is defined in strict analogy to the classical 
action in the deterministic case, but for the limiting procedure that 
needs to be taken in the sense of expectation values, as the sample 
paths of a diffusion process are non differentiable: 



A(t ,ti;q)= / lim At ^ + E 

J to 



2 \AtJ 



dt , (16) 



where E denotes the expectation with respect to the probability den- 
sity p. It can be shown that the mean classical action (16) associated 
to the diffusive kinematics (13) can be cast in the following particu- 
larly appealing Eulerian hydrodynamic form Q: 

P(x,t). (17) 



A{t ,h;q) = ^ dtjd 3 x j (v 2 - u 2 ) - V(x) 



I to 

The stochastic variational principle now follows: the Action is sta- 
tionary, 5A = 0, under smooth variations of the density Sp, and of the 
current velocity 5v, with vanishing boundary conditions at the initial 
and final times, if and only if the current velocity is the gradient of 
some scalar field S(x,t) (the phase): 

mv = VS . (18) 

With the above conditions met, the two coupled nonlinear La- 
grangian equations of motion for the density p (or alternatively for 
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the osmotic velocity it) and for the current velocity v (or alternatively 
for the phase S) are the Hamilton-Jacobi-Madelung equation: 




+ V(x) = 



(19) 



VP 



and the continuity equation: 



dtp 



W[pv\. 



(20) 



By solving equations (19) and (20) the state of the bunch is com- 
pletely determined. Formal linearization of the equations can be 
achieved through the standard De Broglie ansatz yielding the Schroedinger 
equation of the quantum-like models. However, one should bear in 
mind that the real hydrodynamic equations (19)-(20) are the physi- 
cally fundamental objects, while linearizing them to a complex Schroedinger 
equation is a bare mathematical tool that can be useful for calcula- 
tional needs, but bears no physical significance. In particular, the 
complex wave function is devoid of any physical meaning. Thus, in 
the present context, the situation is just the opposite to that in quan- 
tum mechanics, where instead the wave function and the Schroedinger 
equation are the fundamental physical ingredients. 

The observable structure is quite clear: the expectations (first mo- 
ments) of the three components of the current velocity v are the aver- 
age velocities of oscillation of the bunch center along the longitudinal 
and transverse directions. The expectations (first moments) of the 
three components of the process q(t) give the average coordinate of 
the bunch center. The second moments of q{t) allow to determine the 
dispersion (spreading) of the bunch. In the harmonic case, these are 
all the moments that are needed (Gaussian probability density), and 
we have coherent state solutions. In the anharmonic case the coupled 
equations of dynamics may be used to achieve a controlled coherence: 
given a desired state (p, v) the equations of motion (19) and (20) can 
be solved for the external controlling potential V(x, t) that realizes 
the desired state. Lack of space prevents us from commenting further 
on this very important application of our formalism. A thorough and 
detailed study of the controlled coherent evolutions in the framework 
of our stochastic model will be presented in a forthcoming paper ||] . 
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